Reflection from Corners in Rectangular Wave Guides — 
Conformal Transformation* 

By S. O. RICE 

A conformal transformation method is used to obtain approximate expressions 
for the reflection coefficients of sharp corners in rectangular wave guides. The 
transformation carries the bent guide over into a straight guide filled with a non- 
uniform medium. The reflection coefficient of the transformed system can be 
expressed in terms of the solution of an integral equation which may be solved 
approximately by successive substitutions. When the corner angle is small and 
the corner is not truncated the required integrations may be performed and an 
explicit expression obtained for the reflection coefficient. Although applied here 
only to corners, the method has an additional interest in that it is applicable to 
other types of irregularities in rectangular wave guides. 

Introduction 

THE propagation of electromagnetic waves around a rectangular corner 
has been studied in two recent papers, one by Poritsky and Blewett 1 
and the other by Miles 2 . Poritsky and Blewett make use of Schwarz' 
"alternating procedure" in which a sequence of approximations is obtained 
by going back and forth between two overlapping regions. Miles derives 
an equivalent circuit by using solutions of the wave equation in rectangular 
coordinates. Several papers giving experimental results have been pub- 
lished. Of these, we mention one due to Elson 3 who gives values of reflection 
coefficients for various types of corners. 

Here we shall deal with the more general type of corner shown in Fig. 1 
by transforming, conformally, the bent guide (in which the propagation 
"constant" of the dielectric is constant) into a straight guide in which the 
propagation "constant" is a function of position — its greatest deviation 
from the original value being in the vicinity of points corresponding to the 
corner. This type of corner has been chosen for our example because it 
possesses a number of features common to problems which may be treated 
by the transformation method. 

The essentials of the procedure used are due to Routh 4 who studied 
the vibration of a membrane of irregular shape by transforming it into a 
rectangle. After the transformation the density (analogous to the propaga- 
tion constant in the guide) was no longer constant but this disadvantage 
was more than offset by the simplification in shape. 

Until this paper was presented at the Symposium I was unaware of any 

* Presented at the Second Symposium on Applied Mathematics, Cambridge, Mass. 
July 29, 1948. 

1 See list of references at end of paper. 
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other wave guide work based on conformal transformations (as described 
above) except that of Krasnooshkin 6 . At the meeting I learned that 
the transformation method had also been discovered (but not yet published) 
by Levine and by Piloty independently of each other. Levine has studied 
the same corner, see Fig. 1, as is done here. However, his method of 
approach is quite different in that he obtains expressions for the elements 
in the equivalent pi network representing the corner, whereas here the 
reflection coefficient is considered directly. This is discussed in more 
detail at the beginning of Section 6. Piloty's work is closely related to the 
material presented in a companion paper 6 and is discussed in its introduction. 

In this paper the partial differential equation resulting from the trans- 
formation, together with the boundary conditions, is converted into a rather 
complicated integral equation. Numerical work indicates that satisfactory 
values of the reflection coefficient, in which we are primarily interested, 
may be obtained by solving this integral equation by the method of succes- 
sive substitutions. However, the question of convergence is not investigated. 

Although they are here applied only to corners, the equations of Sections 3, 
4 and 5 are quite general. In order to test their generality they were used 
to check the expression for the reflection coefficient of a gentle circular 
bend in a rectangular wave guide, E being in the plane of the bend. The 
work has been omitted because of its length. It was found that the essential 
parts of the transformation may be obtained by regarding the inner and 
outer walls of the guide system as the two plates of a condenser, solving the 
corresponding electrostatic problem (using series of the Fourier type), and 
utilizing the relation between two-dimensional potentials and the theory of 
conformal mapping. 

When the angle of the corner is small we may obtain the series (7-5) 
and (7-11) for the reflection coefficients corresponding to simple (i.e. not 
truncated) E and II corners, respectively (a corner having the electric 
intensity E in the plane of the bend will be called an E corner or an electric 
corner. II corners are defined in a similar manner). When the angle of 
the general E corner shown in Fig. 1 is small we may use the series (7-18). 

The series (7-5) and (7-11) giving the reflection from small angle corners 
are related to the series giving the reflection coefficients for gentle circular 
bends. In fact, if the radii of curvature of the latter be held constant 
while the angle of bend is made small, the series for the circular bends 
reduce to those for the corners. 

As for the limitations of the method, note first that it can be used only 
for wave guide systems in which the dimension normal to the plane of 
transformation is constant throughout. Moreover, the integral equations 
of the present paper, except for the work of Appendix III, are derived 
on the assumption that the dimensions of the guide approach constant 
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values at minus infinity and the same values at plus infinity. When this 
assumption is not met, a conformal transformation may still be used to 
carry the system into a straight guide. However, there appears to be some 
doubt as to the best way of dealing with the resulting partial differential 
equation. One method, discussed in the companion paper , leads to an 
infinite set of ordinary linear differential equations of the second order. 
Again, possibly the Green's functions appearing in Sections 3 and 5 may be 
replaced by suitable approximations. 

1. Representation of Field for Corner or Bend in Rectangular Guide 

Quite often waves in rectangular wave guides are classed as "transverse 
electric" or "transverse magnetic". However, for our purposes it is 
more convenient to class them as "electrically oriented" or "magnetically 
oriented" waves. 8,9 Thus, the electric and magnetic intensities are obtained 
by multiplying 

1 d*A _ dB H = dA 1_ ^B 

iwe dxd{ dy dy iwp dxd£ 

e - L^A + «2 h - - M + J- ** (i-i) 

* um dydt ^ dx ' dx T icctx dyd{ 

Er = - uapA + .— —r Hi = - iweB + — — 

f uotdt 2 uapdr 



by e iut and taking the real part. Here a, /i, and e are the radian frequency, 
the permeability of the medium filling the guide 0* = 1.257 X 10 -6 henries 
per meter for air), and the dielectric constant of the same ( e = 8.854 X 10 -12 
farads per meter for air), respectively. *, y, and f constitute a right-handed 
set of rectangular coordinates in which the f axis is normal to the plane of 
the bend. Equations (1-1) may be verified by substituting them in 
Maxwell's equations. 
The potentials A and B satisfy the wave equation 

d 2 A . d 2 A . d 2 A 2 . 

+ + 7T = O A 

dx' T df T af 2 (1 _ 2) 

a — iw\/ne = i2ir/\{j 

where X is the wave length in free space corresponding to the radian 
frequency oj. 

When the electric vector lies in the plane of the bend, as shown in Fig. 1, 
and the incident wave contains only the dominant mode we set 

.4=0, B = Q sin (irf/a) (1-3) 
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where a is the wide dimension of the rectangular cross-section, the guide walls 
normal to the f axis are at f = and f = a, and Q is a function of * and y 
such that 

S + g-rt.Q-0 

r 10 = ft2*>£*(l - Xoa~V-4)." 2 

The guide walls are assumed to be perfect conductors and hence the tan- 
gential component of E must vanish at the walls. This requires the normal 
derivative of Q to vanish at those walls which are perpendicular to the 
plane of the bend : 

% - 0. (.-5) 

dn 

When the magnetic vector lies in the plane of the bend and the incident 
wave consists of the dominant mode, we set 

A = P, B = (1-6) 

where P is a function of .\ and y such that 

2 2 

^' + f-f - r »» P = °> r »» =- * 2,r Ao (1-7) 

d.v 2 d) 2 

and 

/ J = (1-8) 

at the walls perpendicular to the plane of the bend. In this case the guide 
walls parallel to the plane of the bend are at f = and f = b. 

2. Electric Vector in Plane of Bend 

Figure 1 shows a section of the bend taken parallel to the electric vector. 
b is the narrow dimension of the guide. Let the frequency and the wide 
dimension a of the guide (measured normal to the plane of Fig. 1) be such 
that only the dominant mode is freely propagated. The position of any 
point in this section is specified by the complex number z = x + iy where 
the origin and the orientation of the axes have been chosen somewhat 
arbitrarily. 

The constant k and related propagation constants which appear in the 
formulas dealing with Q and electric bends are given by 

k = (26/Xo) [1 - (X„/2a) 2 ] ,/2 = -irio&A 

7m = m 2 — k 2 ; m = 0, 1, 2, • • • ; 70 = ik (2-1) 

Xo = free space wavelength 
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Since, by assumption, only the dominant mode is freely propagated, k and 
7 m for m > are real and positive. 

We imagine an incident wave of unit amplitude coming down from z B 
in the upper left portion of Fig. 1. What are the amplitudes of the reflected 
wave traveling back toward z 5 and the transmitted wave traveling outward 
to the right towards z 3 ? Our task is to find a Q(x, y), satisfying the wave 
equation (1-4) and the boundary condition (1-5), which represents a 
disturbance of the assumed type. 



z 5 = co e 
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The first step is to find the conformal transformation 

2 = x + iy = f( v + iff) = f(w) (2-2) 

which carries the bent guide (shown in Fig. 1) in the (x, y) plane over into 
the straight guide (shown in Fig. 2) in the (v, 6) plane. This may be done 
by the Schwarz-Christoffel method discussed in Appendix I. This trans- 
formation carries the wave equation (1-4) and the boundary condition 
(1-5) into 

(2-3) 



t9 + tQ + [i + g {v, e)]k 2 Q = o, 



dQ 

dd 



= at 6 = and 6 = tt 



(2-4) 
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where the upper and lower guide walls are carried into = and 6 = ir, 
respectively, and g{v, 0) is given by 

1 + Ah *) = \f(* + M) | 2 T 2 /b z (2-5) 

/ _v _ [ch v + cos 6] 2a 

gKV ' ) '' = [ch(v - t) - cos dY\ch{v + t) - cos 6]° ~ { ' 

Here ch denotes the hyperbolic cosine, f'(v + id) denotes the first deriva- 
tive off (to), and from Appendix I, 2ira is the total angle of the bend. / is a 
parameter which depends upon a and the ratio d/d where d = | z 4 — z 1 
and d = \ z 4 — z 6 | in Fig. 1. A table giving values of / for a 90° bend 
(a = 1/4) appears in Appendix I. 

That the propagation constant is no longer uniform in the transformed 
guide shows up through the fact that the coefficient of k 2 Q in (2-3) is now a 
function of the coordinates (v, d). g(v, 9) measures the deviation of the 
propagation constant from its value at n = — oo . For example, if we 
consider a wave front coming down from z 6 we expect it to get past z 4 before 
it reaches z . In Fig. 2 the same wave front is tilted forward corresponding 
to a high phase-velocity (or small propagation constant) at z 4 where v = 
and 6 = tt. This is in line with the fact that the coefficient of k 2 Q in (2-3) 
vanishes at z i by virtue of (2-6). Similar considerations hold at z x and z 2 . 

What is our reflection problem in terms of the transformed guide? In 
addition to satisfying the two equations (2-3) and (2-4) Q must behave 
properly at infinity. For large negative values of v, Q must represent an 
incident wave plus a reflected wave. The incident wave is of unit amplitude 
and the reflected wave is of the, as yet, unknown value R E . For large 
positive values of v Q must represent an outgoing wave. Thus Q must 
also satisfy the two equations 

Q = e -ikv _|_ R se ikv t v _> _ oo (2-7) 

Q = T B e~ ik * , v -> « (2-8) 

where the subscript E appears on the "reflection coefficient" R E and the 
"transmission coefficient" T s to indicate that here we are dealing with an 
electric corner. 

Our problem is now to take the four equations (2-3, 4, 7, 8) and somehow 
or other obtain the value of R E . We are not so much interested in T E 
because it does not have the practical importance of the reflection coefficient. 
There are at least two different ways we may proceed from here. One 
is to transform the differential equation plus the boundary conditions 
into an integral equation which may be solved approximately by iteration. 
Another way is to assume Q to be a Fourier cosine series in 6 whose co- 
efficients are functions of v. Substitution of the assumed series in the 
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differential equation (2-3) gives rise to a set of ordinary differential equations 
having v as the independent variable and the coefficients as the dependent 
variables. The integral equation method is used in this paper. The 
second method is discussed in the companion paper. 

3. Conversion of Differential Equation into an Integral Equation 

The differential equation (2-3) may be converted into an integral equation 

by using the appropriate Green's function in the conventional manner. 

The only modifications necessary are essentially those given by Poritsky 

and Blewett 1 in a similar procedure. 

The conversion is based upon Green's theorem in the form 

where the integration on the rightextends over the rectangular regional < v<v 2 , 
< < 7T (inside the straight guide associated with (v, 0), i.e. the guide of 
Fig. 2) except for a very small circle surrounding the point (v , B ). 
G = G(v , 0o ; v, 0) is the Green's function corresponding to 

d ll + *X + k 2 V = (3-2) 

dv* T 60 2 

in the region — * < v < «. , < 0<7r subject to the boundary condition 
dV/dn = on the walls (dV/dd = at = and = *■). G becomes 
infinite as —log r when r — > 0, r being the distance between the variable 
point (», 0) and the fixed point (v , O ). Poritsky and Blewett* have shown 
that, in the notation (2-1), 

G = E e m7 Z l cos md cos mef^^ (3-3) 

m=-0 

e = 1, e m = 2 for m = 1, 2, 3 • • • 
Equation (3-1) leads to 

= k 2 r dv [ dd g(v, d)QG 

Jvi Jo 

from which the required integral equation for Q is found to be 

Q(v , 0o) - 

,-*•• + ^ / '"dv fdd g(v, 6)Q(v, 0) it*»7~ c ° s m0 o cos md r 1 "- ' 1 '" (3-5) 
2irJ—co Jo "»=o 

* We have replaced their i by —i since here we assume the time to enter through the 
factor e iut instead of e~ iut . 



(3-4) 
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where y m is given by (2-1). The term e~ ihv o comes from the first integral 
on the left side of (3-4) as i'i — > — oo . Equation (3-5) is a general equation 
which may be applied to a number of wave guide problems by choosing a 
suitable function g(v, 6). For the corner of Fig. 1 g(v, 6) is given by (2-6). 
If g(v, 0) approaches zero when | v | becomes large, as it does for the 
corner, expressions for the reflection coefficient R E and the amplitude T E 
of the transmitted wave may be obtained by letting v — > ±°° in (3-5). 
For very large values of | v | the contributions of all the terms in the summa- 
tion except the first (m = 0) vanish. Comparison of the resulting expression 
for Q(v , O ) with the limiting forms (2-7) and (2-8) defining R B and T B gives 

Rm = - £ £ dv £ dd g(v, 6) Q(i, 9)e- ikv (3-6) 

Tm = 1 - £ £ dv j[" dd g(v, d) Q(v, e)e tk ° (3-7) 

Since the integrands involve the as yet unknown Q(v, 6) these expressions 
are not immediately applicable. In fact, if we knew Q(v, 6) it would not be 
necessary to use these integrals for R E and T E — we could simply let v — » ± oo 
and use (2-7) and (2-8). Nevertheless, (3-6) and (3-7) are useful in obtain- 
ing approximations to R E and T E when approximations to Q are known. 

In Appendix IV it is shown that R E is the stationary value, with respect to 
variations of the function Q, of an expression made up of integrals containing 
Q in their integrands. From the integral equation it follows that when 
k — ► 0, i.e., when the frequency decreases toward the cut-off frequency of the 
dominant mode, Q becomes approximately exp (— ikv). Furthermore, R E 
approaches zero. This is in contrast to the apparent behavior of R H which, 
according to the discussion given in Section 5, may possibly approach — 1 
under the same circumstances. Thus reflections from the two types of 
corners, or more generally, irregularities in the E plane and in the H plane, 
appear to behave quite differently as the cut-off frequency is approached. 

R E and T E are not independent. Since the energy in the incident wave is 
equal to the sum of the energies in the reflected and transmitted waves we 
expect 

ReR*b + T E T* e = 1, (3-8) 

where the asterisk, denotes the conjugate complex quantity. In addition, 
there is a relation between R E and T E which for a symmetrical irregularity, 
i.e. for g(v, 6) an even function of v, states that the phase of R E differs from 
that T B by ±7r/2. In this special case T B is determined to within a plus or 
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minus sign when R E is given. These relations may be proved by substituting 
various solutions of equation (2-3) for Q and Q in the equation 

\q£-Q»£\ =\ Q d 4-Q d -§] (3-9) 

where m and v 2 are large enough (»j negative and v 2 positive) to ensure that 
Q and (5 have reduced to exponential functions of v. Equation (3-9) follows 
from Green's theorem. When Q is taken to be the solution for which (2-7) 
and (2-8) holds and Q its conjugate complex Q*, equation (3-8) is obtained. 
Keeping the same solution for Q but now letting Q denote the solution 
corresponding to an incident wave of unit amplitude coming in from the 
right: 

Q l = e ikv + R 1 g-»' ) v — > oo 

Qi = T!e ikv , »-*-*> 

gives T = Ti where we have dropped the subscript E and have assumed 
that g(v, 0) may be unsymmetrical. Taking Q to be Q x gives 

Rf$ + RtT = 

which is the relation sought. In the symmetrical case R = Ri, R/T + R*/T* 
is zero and hence R/T is purely imaginary as was mentioned above. The 
same relations hold for R B and T B . These results are special cases of a 
more general result which states that the "scattering matrix" is symmetrical 
and unitary for a lossless junction. 10 

4. Approximate Solution of Integral Equation 

A first approximation to the solution of the integral equation (3-5) is 
obtained when we assume that the non-uniformity of the propagation 
constant has no effect on Q. Thus we put 

Q(»( V} e) = e~ ikv (4-1) 

in the integral on the right and obtain an expression for the second approxi- 
mation Q l2) (v, 0)> and so on - Here we sha11 not S° bevond Q {2) ( v > e )- 
It is convenient to expand g(v, 6) in a Fourier cosine series 

CO 

g(v, 6) = £ a n (v) cos nd 

(4-2) 
a n (v) = - [ g(v, 6) cos nd dd, e = 1 ; €„ = 2, n > 0. 

IT Jo 
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The second approximation, obtained by substituting (4-1) in (3-5), may then 
be written as 

00 

Q (2) (v , do) = e~ ikv ° 4- k 2 2~ l E 7m 1 cos mdo 

(4-3) 
a m (v)e- ikv ~ ]v - Vohm dv. 



L 



The nth approximation R E n to the reflection coefficient (when the 
electric vector lies in the plane of the bend) is defined in terms of Q in) by 



Limit Q {n \v, d) = «-'*' + lft"V" (4-4) 

V-* — CO 

R B n) is also equal to the integral obtained by replacing Q in (3-6) by Q (n-1) . 
We have 

rP = 0, R B 2) = -ik2~ l r ao(v)e- 2ikv dv, 

J— to 

gp = R B > ) -ik 3 E(4 7m e m )- 1 (4-5) 

m— 

• [ dv a m (v ) f dva m (v)e- ik(v+Vo) - 

J- oo •*— to 



)-lt.-» |- 



where y m is given by (2-1). 

The results of this section have the same generality as the integral 
equation (3-5) in that they are not restricted to corners. 

5. Truncated Corner — Magnetic Vector in Plane of Bend 

When the magnetic vector lies in the plane of the bend the reflection 
may be calculated by a similar procedure. The wide dimension a of the 
wave guide now replaces the narrow dimension b in Fig. 1. We shall call 
the result of making this change the "modified Fig. 1". We again assume 
the frequency to be such that only the dominant mode is propagated without 
attenuation. In place of equations (1-3, 4, 5) involving Q we have those of 
(1-6, 7, 8) involving P. 

The conformal transformation which carries the modified Fig. 1 into 
Fig. 2 leads to 



d'P . d'P . t . , , -M 2 



+ ±L + [i + g ( V) e)h 2 p = o 

dv 2 dd z 8 ' (54) 

P = at 6 = and 6 = t 
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where 

k = 2a/X - -flWx, c = (k 2 - 1) 1/J = ak/b 

5 2 m = m* - 1 - c 2 = w 2 - * 2 , 81 - iV (5-2) 

Xo = free space wave length, m = 1, 2, 3 ■ ■ • 

and 

|/Ld(* + id) \***/a* = 1 + g(v, 6). (5-3) 

Here f moA (w) pertains to the modified Fig. 1. Since the expression for 
f'(w) given in Appendix I is proportional to b and since the modified trans- 
formation contains a in place of b, it follows that g(v, 6) for the magnetic 
corner is exactly the same function, given by (2-6), as for the electric corner. 
It is again assumed that the incident wave coming down from the left 
in the modified Fig. 1 is of unit amplitude and of the dominant mode. 
At large distances from the corner 

P = [e~ icv + R„ e icv ] sin 6, v -» - 00 

(5-4) 

p = T„ e~ ,cv sin 6, v -> + =c 

which serve to define the coefficients of reflection and transmission. The 
subscript H on the reflection and transmission coefficients indicate that 
here we are dealing with a magnetic corner. 

The conversion of the differential equation into the integral equation now 
employs the Green's function 

G = 2 f) 8~ l sin m.6 sin m0<T U " "" (5-5) 

m=l 

which corresponds to 

q + q + tv-o 

dv 2 90 2 
V = at 6 = and d = ir 

The integral equation for P is found to be 

P(v , 0o) = e~ icv ° sin do 

+ — [ dv I dd g(v,e)P(v,$) '£,28^sinm0 smm6e~ u " ol 
2tt J- * Jo m-l 

where the parameters are given by (5-2). This is a general equation. 
For the corner of the modified Fig. 1 g(v, 0) is given by (2-6). 



(5-6) 
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By letting v — > — * we obtain the exact expression 

• 2 * oo »x 

Rh = -— dv / </0 gfo 0)P(j., 0)6-'" sin (5-7) 

ttc J-x Jo 

When dealing with the electric corner we saw that R E — > as k — » 0. 
The presence of c in the denominator of (5-7) suggests the possibility that 
Rh — > — 1 as c — * 0. For P w must remain finite and this may perhaps come 
about through P(v, 8) — * in the region, say around v = 0, where g(v, 6) is 
appreciably different from zero. This and the fact that P(v, d) must 
contain a unit incident wave suggest that for v < the dominant portion of 
P(y, 6) is 2/ sin cv which gives P w = — 1. Incidentally, it is apparent that 
the approximations for P(v, 6) given below in (5-8) and (5-10) (and therefore 
also the approximations (5-11) for R H ) fail when c becomes small. 

The first approximation to the solution of the integral equation (5-6) is 

P<«(», 9) = e~ icv sin (5-8) 

When we introduce the coefficients 

2 [" 
bn(v) = - / g(v, 6) sin 6 sin nd dd 

7T JQ 

sin 8g(v, d) = ^ b n (v) sin nd 

n-1 

fo(v) = a (v) — a 2 (v)/2, b n (v) = [a„_i(fl) — a n+ \(v)]/2, n > 1 
we find that the second approximation is 

oo 

P i2) (r , do) = e~ icv ° sin 9 + k2~ 1 £ 5" 1 sin m$ 

• f *.(») 

J— 00 

The successive approximations to the reflection coefficient are 

Ki, l) = 0, P< 2) = -f f +C °dvbMe- 2i " 

2c J- * 

itf* = Pj, 2) - u- 4 £ (4CU" 1 [^ dv b m (vo) (5-11) 



(5-9) 



(5-10) 



k _ lc „_i„_ l , u m 



£ 



+ 00 

— ic(p+v )— I •■ — » o 1 5 r 



110 BELL SYSTEM TECHNICAL JOURNAL 

6. Series for R {2) When Corner Has No Truncation 

The integrals which appear in the approximations for the reflection 
coefficients are difficult to evaluate in general. This section serves to put 
on record several expressions which have been obtained for i? (2) when the 
corner is not truncated. Corresponding evaluations of R w would be 
welcome since the work of Section 7 for small angle corners indicates that 
2J(3) _ _R(2) i s f the same order as R m . However, I have been unable to 
go much beyond the results shown here. 

As mentioned in the introduction, H. Levine has studied the effect of a 
corner in a wave guide by representing it as an equivalent pi network having 
an inductance for the series element and two equal condensers for the shunt 
elements. Early in 1947 he derived the following expressions (in our 
notation) for the elements corresponding to a simpe E corner:* 

vr.-*[*(nn) -*<-*>] 

B b /Y = (for) -1 cot Gftr/2) 

where Y is the characteristic admittance of the straight guide, iB a the 
admittance of one of the two equal shunt condensers, -iB b the admittance 
of the series inductance, ^(x) the logarithmic derivative of T(x + 1), and 
fiw is the total angle of the simple corner (for no truncation we set = 2a). 

When the reflection coefficient for the corner is computed from the 
equivalent network for the case — > it is found to lie between the approxi- 
mate value R E 2) given by (7-3) and the considerably more accurate value 
R B Z) given by (7-5). All three approximations are of the form A(? + 0(/3 3 ) 
where A differs from approximation to approximation but is independent of 
/3, and 0(/3 3 ) denotes correction terms of order /3 3 . Since R E 3) gives the exact 
value of A, it may be regarded as the standard when the three approxima- 
tions are compared. If this comparison be taken as a guide, it suggests 
that the rather cumbersome expressions (6-2) and (6-5) for Rb given below 
are not as accurate as the simpler expressions resulting from Levine's work. 
Dr. Levine has also obtained corresponding results for the general £-corner 
of Fig. 1. It is hoped that his work will be published soon. 

When the corner is not truncated it is convenient, as mentioned above, 
to replace 2a by /3 so that /3tt is the total angle of the bend. For no trunca- 
tion / = and (2-6) becomes 



g(v, 8) = 



chv -f- cos 6 
chv — cos 6 



- 1. (6-1) 



* I am indebted to Dr. Levine for communicating these expressions to me. 
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From (4-2) and (4-5), or from (3-6), 

*P- -~ I' dv f d8g(v,0)e-" 

ZlT J— oo Jo 

= - ik f T (n - ik)T{ n + ik) v (- 2 fl»-(ffl-- 

£i n !(« - 1 ) ! 2 £$> (2m) ! (» - ».) ! 

where we have expanded g(t>, 5) as given by (6-1) in powers of 
cos 6/ch v and integrated termwise. The notation is (a) — 1, (a) n = 
a(a + 1) • • • (a + n - 1). 

For a right angle corner /3 = 1/2, and a more rapidly convergent series 
may be obtained by subtracting the sum of the series corresponding to 
k = 0, namely 



log 2 = E ^ (6-3) 

n -i n\2n v ' 



(1/2), 

n— 1 

Thus for /3 = 1/2 

n-1 " (6-4) 

^1 = n-A/sinh **, -4n = /li IT (1 + k 2 m~ 2 ), n > 1 

The rate of convergence of the more general series (6-2) may be increased 
in a somewhat similar way. It is found that 

Bp - -j[j - 2^(1 - Ad - £ (2 + «(1 - A,) 

- j| (23 + 2O0 2 + 20 4 )( 1 -.43) J 

7 ~ K + Z, (6-5) 



J (2*»)! £2, (« - m)\n H £?, (l/2) m m(m - 0) 

where .5772 • • • is Euler's constant, ty(x) is the logarithmic derivative of 
Il(x) = T(x + 1), and A n is given by (6-4). 
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The results corresponding to R ( ,P are quite similar. When the corner 
is not truncated 

(2 ) _ • 2 v r( rc - ic)T(n + ic) y> (-2j8) 2 mG3)»-» 



£i (« + 1)1 (« • - l)12c£o (2m)\(n - m)\ 

, _*£ [j - /3 2 (1 - i,) - £ (2 + /3 2 )(1 - A 2 ) 
2c |_ ^ 



(6-6) 



-^(23 + 2O/3 2 -f-20 4 )(l -i 8 ) - 



/ = 1 - ¥(1 - 0) - .5772 

,. . f, (1/2 ~ 0)m 

- 0U - P) Z, ( 1/2 ) mW ( w _ ^(m - + 1) 

in which A n is obtained by replacing c by ft in the expression (6-4) for A n . 
The evaluation of the integrals for R^ and R% } for general values of / 
appears to be difficult although it is possible to obtain approximate expres- 
sions for the case when / is large. 

7. Reflection from Small Angle Corners 

The expressions for R {2) and R m may be evaluated approximately when 
the angle of the corner is small. It turns out that, for t = 0, they are of the 
same order of magnitude and both of them must be considered. Moreover 
R M for n > 3 differs from R (z) by terms of the same order as those neglected 
in our approximations so that there is no point in going to the higher values 

of n. 

We first obtain the approximation for R E for a corner with no truncation 
having the total angle vfi. Since is very small (6-1) may be written as 

g (v, d) = exp \fr] - 1 - fo + 0V/2 ! + 0(/5 3 ) 

tp = log (chv + cos 6) — log (chv — cos 0) 

where 0(0') denotes terms of order /3 3 . The expression <p becomes very large 
near the two points (0, 0) and (0, *■) (the coordinates being (v, 8)). The 
following considerations indicate that this does not invalidate our procedure. 
The remainder, denoted by 0(/3 3 ), in (7-1) is less than | fo | 3 exp | ftp |. 
Near (0, 0) <p is approximately equal to 21og(2/r) where r"- = v°- + 6\ 
Consequently the remainder is less than (2/3 log 2/r) 3 (2/r) 2/J . When the 
expression (7-1) for g(v, d) is set in the integral equation it is seen that all 
terms, and in particular the remainder term (by virtue of the inequality 
just stated), of the double integral converge at (0, 0). Hence the contribu- 
tion of the remainder term is of order (3 3 even in the worst case when the 
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Green's function is replaced by —log r. A similar result holds for the other 
point in question, namely (0, ir). 

Integrating (7-1) from 6 = to = t and using equations (A2-1, 3) of 
Appendix II gives 

(h(v) = P[h(v, v) - I,{v, v)] +■ O03 3 ) 

= 4/3 2 £ »~* e~ 2nv + O03 3 ) (7-2) 

n— 1,9,5.* • • 

where v > 0. We consider only positive values of v since g(v, 0) and the 
a n (v)'s are even functions of v. Thus (4-5) yields 

R? = -ik2(? 2 » _1 (» 2 + * 2 ) -1 + 0(/3 3 ) (7-3) 

n— 1,8,8, ••■ 

This is an approximation to the exact value given by the double series in 
(6-2). Comparison of (6-5) and (7-3) when /3 and k approach zero gives, 
incidentally, 



E»- 2 t(»- 1/2)- 1 = H,m- 1 



From (4-2), (7-1) and the expansions (A2-2) of log (chv ± cos 0) it follows 
that 

a m (v) = 4pnr 1 <r m \'\ + 0(/3 2 ), m = 1, 3, 5, • ■ • 

(7-4) 
a m (v) = 0(/3 2 ) , m = 0, 2, 4, 6 • • • 

Equations (4-5), (A2-4), the relation y 2 m = w 2 — & 2 , and (A2-8) give us the 
answer we seek : 

R { P = R^ - ik*2& E T» m~ 2 J(m, m, k, y m ,0,0) + 0(/3 3 ) 

m=l,3,5.--- 

= -iklf? £ 7 « l «"* + 0C8 8 ) (7-5) 

m=l,3,5.- • ■ 

It is not necessary to go to R B because it differs from Rg* by only 0(/3 3 ). 

When H lies in the plane of the bend the reflection from a small angle 
corner with no truncation may be obtained by much the same procedure. 
For brevity we shall not write down the order of magnitude of the remainder 
terms. From (5-9), (A2-1), and (A2-3) 

h(v) = fl„(r) - a 2 (v)/2 (7-6) 

= /r°[/, - h- (7 3 - Id/2) 

= /r[2e- 2 " + 4 £ ,rV 2n " - 4 Z (n - D"*«"*l 
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where we have written /„, for I m (v, v) and assumed v > 0. Then, using 
(5-11), 

sg> = -itc*ti** + 2 e »~v + *r* 

•-»•*•- (7-7) 

-2 Z »(» t -i)"V + ^)" 1 ] 

n-2.4.-" 

When we put 
6,(1.) = [fln-i(») - On+iWl/2, n > 1 (7-8) 

&„(*) = 2/3[(« - 1)-V«-*>W - (« + l)- 1 e- (n+1) ' , '>], » = 2, 4, 6, • • • 

K(v) = 0(/3 2 ), » = 1, 3, 5, • • • 
in (5-11) and use the results of Appendix II we obtain 

n-2.4.6--- 

+ (n + l)~ 2 -/(» + 1, » + 1, c, 5„ , 0, 0) (7-9) 

-2(n 2 - ir 1 7(w-l,«+ l,c,8 B ,0,0)] 

The values of the first two J's, obtained by setting m = n ± 1 in (A2-7), 
may be simplified by using 

c 2 + (» ± 1 + 5) 2 = 2(» ± 1)(» + 5) 

where we have dropped the subcript n from 5„ . In order to eliminate 5 
from the denominator we multiply both numerator and denominator by 
n — 5 and use 

(n - 8) (8 + In ± 2) = (» ± l) 2 + c 2 - 8(w ± 2) 
M 2 _ 5 2 = i + C 2 = K 2 

Setting in the value, given by (A2-9), of the last / and separating the 
terms (into those which contain the first power of 5 and those which do not) 
enable us to write the term within the square brackets in (7-9) as 

4« 2 5„T (» - 1) - 1 , ( n + 1) + 1 

**(»*- 1)« k 2 L(» - dv + (« - 1) 2 } (» + ov + (» + 1) 2 } 

, 2n{2(n 2 -r-c 2 )-KV-l)} ] f? Q) 

"*" K 2 (» 2 - 1) 2 (« 2 + C 2 ) J 

It is found that when (7-10) is put in (7-9), the contribution of the first 
two terms within the square bracket of (7-10) exactly cancels the summation 
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which is taken over 3, 5, 7 • • • in the expression (7-7) for Rh\ Moreover, 
if we make use of 

£ Mn - I)" 2 = 1 

n-2,4.6,--- 

we see that the contribution of the last term within the square brackets of 
(7-10) cancels the remaining terms in Ri?\ Only the contribution of the 
first term in (7-10) remains and it gives 

RP = -itfc- 1 £ n 2 (n 2 - 1)- 2 C + 0(/3 3 ) (7-11) 

n=2.4,8.- •• 

The relative simplicity of this result indicates that there may be another 
method of derivation which avoids the lengthy algebra of our method. 

Recently approximate expressions for the reflection coefficient of gentle 
circular bends have been published 7 . In our present notation these may be 
written as 



x. - - avfe - «* L 

|_ /4 m-l.3.6 



cos u - e-*"""" 1 ' 



ro-1.3.6.-- TflWYn 



Rg m - ia 2 Pl 2 



["sin u _ y cos « - g ulJe n 2 "1 

L87T 2 C 2 ' „- 2 il.... v*Cdn (« 2 - 1) 3 J 



where /37r is the angle of the bend, pi is the radius of curvature of the center 
line of the guide and u is 2v times the length of the center line in the bend 
divided by the wavelength in the guide: 

u = fiirkpi/b = (3ir 2 cpi/a 

The first expression for u is to be used in R B and the second in R H . If we 
now let /3 — > 0, keeping pi fixed, then u — > 0. The trigonometric and expo- 
nential terms may be approximated by the first few terms in their power 
series expansions, and part of the series which make their appearance may be 
replaced by their sums given, for example, by equations (4.1-7) and (4. 1-8) 
of reference 7 . After some cancellation, the above expression for R B and R n , 
which hold for gentle circular bends, reduce to (7-5) and (7-11), respectively, 
which hold for the sharp corners. In other words, the reflection coefficients 
for both the sharp and the circular bends approach zero as /3 — > 0, and 
furthermore their ratio approaches unity. 

We shall merely outline the derivation of the approximation R B 3) for a 
truncated corner. Instead of (7-1) we have from (2-6), 

g(v, 0) = exp M - 1 - cup + <*V 2 /2! + 0(<* 3 ), (7-12) 

ip = 2 \og[chv + cos 6] — log [ch(v — t) — cos 6] — log [ch(v -f t) — cos 6] 
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The Fourier coefficients of g(v, 6) may be obtained by using the results of 
Appendix II. Assuming v > 0, m > 0, 

a (v) = 2a(v - l)f(v) + 2~W{Ah{v, v) + h(v - t, v - t) 

+ h(v + t, v + t) - 4h(v - t, v) - 4/ 2 (f +/,») + 2 /i(» - /, » + /) }, 
a m (») = 2aw- I [-2(-) m e- m l"l + e — 1— *l + e -"l»+t|] (7-13) 

where ^(») = 1 when < v < t and ^(w) = when v > t. Substitution 
of the values (A2-3) for h and h gives 

0o(v) = [2a(D - + 2a 2 (» - 0WW C 7 " 14 ) 

+ a 2 Z «- 2 [4e- 2 "" + e - 2 ""- 2 "' - 4(-)«e- 2n "- n ' 

n— 1 

_^_ g-2n|p-«| _j_ 2«-"l' -l l-*l»+'l — 4(— ) n e- n \ v ~ t \- nv ] 
The second approximation to the reflection coefficient is 
rW = jafr-i s { n 2 kl _ ia*k- 2 2- l (2kt - sin 2kt) 

- ika*J2 »-V + k 2 )~ l {2-(-) n 2e-< (7-15) 

n=l 

+ [1 - 2(-) n e-' u + <r 2n ']cos 2kt 
+ wr 1 ^" 2 "' - (-)»2e-"']sin 2*/} 
The typical term in the summation (4-5) for R B is 

_ * f + " ^ fc W f A a m (,)«- ,A ' ( " +ro) - |r - rol7m (7-16) 

When w = 0, eo = i, 7o - »*, and a (v) is 2a(u - + 0(a 2 ) for < v < t 
and is 0(a 2 ) for V > /. The integral may then be approximated by replacing 
the upper limit °o in (A2-14) by /. The value of (7-16) for m = is found 
to be, to within 0(a 2 ), 

2-W(<r"*' - 1) - (3/4)*a**-*(sin 2kt - 2kt) (7-17) 

When m > 0, e m = 2, yl = w 2 - #, and the substitution of the value 
(7-13) for a m (v) enables us to express (7-16) as the sum of six 7's where / is 
defined by (A2-4). The J's may be evaluated with the help of (A2-7) and 
(A2-8). Substitution of this value of (7-16) and the value (7-17) for m = 0, 
together with R^ given by (7-15), in the expression (4-5) for R^ gives 
our final result 

R^ = ialThm* kt + a 2 / 2 2- 1 (e- 2ifcl - 1) (7-18) 

CO 

+ io? [4" 1 ife- 2 (2/fe/ - sin 2kt) - -Bsin 2kt + k £ »"*y?4j, 
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where 

00 

B = Y. n-' l [e- 2nt - 2(- )««-"] 

n = l 

An = cos 2kt - [2cos kt - (-)V 7 "'] 2 

Equation (7-18) is an approximation, to within terms of order a 2 , for the 
reflection coefficient of a truncated corner which turns through a small 
angle 2mc. The electric vector lies in the plane of the bend. When t = 0, 
(7-18) reduces to (7-5) by virtue of 2a = 0. 

APPENDIX I 

CONFORMAL TRANSFORMATION OF TRUNCATED CORNER 

We shall use a Schwarz-Christoffel transformation* to carry the guide of 
Fig. 1 into the straight guide of Fig. 2. The first step is to transform the 
interior of Fig. 1 into the upper half of an auxiliary complex plane which we 
shall denote by f. Let the points 2i , z 2 , z 3 , z 4 , z 5 in Fig. 1 correspond to 
the points - h, h, 1, », -1 in the £ plane. A suitable transformation 
is then 

z = D + E [\ T . + ft)-( T - A)-( T - l)-i(r + l)-yr (Al-1) 
Jo 

where D, E and h are to be determined from the geometry of Fig. 1. Because 
of the symmetry of our transformation about the line joining z and z 4 
it follows that z = z corresponds to f = 0. Hence D = z . As f travels 
from 1 — e to 1 + e, e being very small and positive, along a semicircular 
indentation above f = 1, z as given by (Al-1) increases by 



|2\-o 



£(1 - fl-r 1 ^ f (r - I)" ir = =^(1 - ^) 

while, according to Fig. 1, it increases from » + t'O to » + i6. Hence we 
set the real part of E equal to -26ir" 1 (l - #)"• We have tacitly assumed 
the factors in (Al-1) to have their principal values at t = 1 + t and also 
that < h < 1. As z goes from Si to z 2 , f goes from -// to +h. In this 
range arg(r + h) =» and arg(r - fr) = t. 
Consequently, if | z 2 — Zi | = f, then 



z 2 — Zi = fo" 



= -Ee~ ia 'l (If- - ^""(l - r 2 )- 1 ^ 

J-/i 



* See, for example, S. A. Schelkunoff, Electromagnetic Waves, New York (1943) 
pp. 184-^187. 
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and we see that E is purely real. Hence 

f = 2far- 1 (l - W)« ( (# - r 2 ) _0 (l - t 2 )- 1 dr 

J-h 

is an equation from which h may be determined as a function of I. Setting 
t 2 = /f 2 .v, expanding (1 — h 2 x)~ l in powers of h 2 and integrating termwise 
leads to 

L = T(T^r (1 " w***^ > ;l - a]lii) 

= T ra-"«) g) //1_2ajF( " - * * -«;*-«; « (ai-2) 



sin 7ra T(f — a) 

where we have used relations from the theory of hypergeometric functions. 
The term 1/sin wa is the reduced form of an original term containing a 
hypergeometric function which has been evaluated by the binomial theorem. 
The second and third expressions are suited to calculation when h 2 < 1/2 
and h 2 > 1/2 respectively. 

Now that the guide of Fig. 1 has been transformed into the upper half of 
the £ plane, the next step is to transform this upper half into the straight 
guide of Fig. 2. We want £ = — 1, i.e. Z5 , to go into v = — 00 and f = 1, 
i.e. 23 , to go into v =■ + 00 . Again using the Schwarz-Christoffel formula 
with w = v + id (the exterior angles at v = ± <*> are equal to ir) 



w 



- Dx + Er f (r + lr'G- - I)' 1 dr (Al-3) 

Jo 



We take the point Zq in Fig. 1 to correspond to v = 0, 8 = in Fig. 2. Since 
this corresponds to f = 0, Di must be zero. Also dw/d£ is real because w 
traverses the walls of the guide of Fig. 2 as f moves along the real axis in the 
f plane. Hence Ei is real. As f goes from 1 — e to 1 -f- e around a small 
circular indentation above f = 1, to changes from co to 00 -f ix. Thus 

iTT = Ei^-wr) or Ex = -2 (A 1-4) 

When (Al-3) is integrated, (Al-4) inserted, and the result solved for f 
we obtain 

£ = tanh to/2 (Al-5) 



• 
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The function we require is obtained by differentiating (Al-1) and (Al-3): 

= £(f - rt-tf - lr^cr 2 - i) 
= a - h 2 n? - 1?)-%/* 

= bf ch 2 w/2 T 

jr lshi(w - t) sh\(w + OJ 

b[ (e" + l) a _T 

DJ 



(Al-6) 



w M*- 1 - l)(e u,+< - 



where 



h = tanh 1/2 (Al-7) 

For a 90 degree corner a = 1/4 and 

i = 2 1/2 (1 - rf/do) (Al-8) 

where, in Fig. 1, d = \ Z\ — z 1 and do = \ Zi — Ze |. In order to obtain 
the relation between t, defined by (Al-7), and d/do various values of h? 
were picked and the corresponding values of / and d/do (using (A 1-2) and 
(Al-8)) computed. Representative values are given in the following table. 



d/do 


/ 


d/do 


t 


1.000 





.5796 


1.2302 


.9041 


.0633 


.5385 


1.4910 


.8565 


.1417 


.5000 


1.7594 


.8292 


.2007 


.4615 


2.0634 


.7745 


.3500 


.3727 


2.8872 


.7196 


.5421 


.2804 


4.0096 


.6919 


.6549 


.1708 


5.987 


.6273 


.9624 


.0959 


8.294 



APPENDIX II 

Integrals Associated with Corners of Small Angle 

The derivation of the integrals encountered in Sections 7 and 8 will be 
outlined here. The first ones are 

Ii(u, v) = I los(c//. n — cos 6) log (ch v — cos 0) dd 

IT JO 
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h(u, v) = - / logfcA W - cos 0) I03 {ch v + cos 5) </0 

* J ° (A2-1) 

/,(« t .) = _ / cos 20 log (ch u - cos 0) log (ch v - cos 0) dO 

IT JO 

U(u, v) = - I cos 20 log (ch u - cos 0) log (cA r + cos 0) dd 

TT JO 

Assuming u and v to be positive and using the expansions 

CO 

log(r// h - cos 0) = log(e"/2) - 2^ /r'^'^cos w0 

(A2-2, 

log(f// u + cos 0) = log(e u /2) - 22 (-) n «" l «- , "cos nB 
leads to 

/,(«, ») = log(««/2) log(e"/2) + 2E n-*-e-"<-"- 

n— 1 

CO 

/,(«, ») = log(e-/2) log(«»/2) + 2Z (-)irV-™ 

n=l 

/ 3 («, ») = -<r 2u log(ey2) - e - 2v \og(e"/2) + 2<r«- 

(A2-3) 
+ 2£ *-»(* + 2)- 1 e-""-" , '(<r 2u + <r 2 ") 



n=l 



/ 4 (m, v) = -<rMog(e"/2) - «- 2 "log (e u /2) - 2«--« 

CO 

+ 2£ (-)■*-»(« + 2)- 1 e-»"-» (<r a " + «-»•) 

When « or v are negative they are to be replaced by their absolute values 
in the expressions (A2-2, 3). 

Now we consider the double integral 

• dT °L dV (A2-4) 

•exp [ — n 1 7' — r I — m \v — s\ —ic(v + v ) — 5 | v — i'o | ] 

in which n, m, c, 8 are real and positive and r and s are real. The double 
integral may be reduced to a single integral by substituting 

e 41 — •« = * f + " (5 2 + rr 2 )-y*<"-"o> rf:r> ( A2 -5) 

7T J-w 
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interchanging the order of integration, and integrating with respect to v 
and vo . Assuming s — r > 0, the integral is then evaluated by closing the 
path of integration by an infinite semicircle in the upper half plane and 
calculating the residues of the integrand at the poles id, c -\- im, — c -f- in'. 



Jin, m,c,8; r, s) = / 

J-oe IT 



idfinie ax 



- tt(6 2 + .r 2 )[ M 2 + (x + c) 2 ][w 2 + (x - c) 2 ] 

[-(6+ic)i+(t-ic)r 

,-^-^r (A2-6) 

+ 



m[8°- + (c + m») 2 ][m 2 + (2c + im) 2 } 

111— 0»+2»c)» 



+ , 



" <» 2 ]J 



M [5 2 + (c - / M ) 2 ][w 2 + (2c 

Substituting special values for the parameters gives the results required 
in the text. Thus, 

J(m, m, k,y; I, I) = c~ 2ikl J(tn, m, k, y; 0, 0) 

J(m, m, k, y; — /, /) = e 2ikl J(m, m, k, 7; 0, 2/) 

J(m, m, k,y;-t,0) = e 2ikt /(*», m, k, 7; 0, /) (A2-7) 

7(w, m, c, 5; 0, 0) = > 2 , nr< . 2 , , , .v 2 . 
(c 2 + m 2 )[c- + (w + 5)-J 

which hold irrespective of any relations between the parameters. The 
derivation of the last result is simplified by setting a. = c -\- im, a = c — im 
and factoring the denominators in (A2-6) so as to obtain terms of the 
form a ± iS, a ± iS. 

When 7 2 = m 2 — k' 2 considerable simplification is possible and we obtain 



/ („, M , i)7; o,o)-xg-- r -»-J 

,, , n A ye~ l ' |~e~ 7 ' e~""(m cos kl — k sin &/) 
7K m, k; y, 0, /) = - \-~ ^T+^- 



(A2-8) 



If we put ft = n — 1, m = ** + 1, and set 5 2 = »* — 1 — c 2 = /; 2 — k 2 

where k 2 = 1 + c 2 , (A2-6) yields, after some reduction, 

./(» - 1, » + 1, c, 6; 0, 0) = 7 — --- + 



{en + iS)* 2(» + 1)(1 - ic) 2 (» - ic) 

± P?__ 

- ic) 2 (n + ?'c) 

»5[2(« 2 + c 1 ) - *V - 1)] 



+ ( "_+- 1)6 - (A2-9) 

T 2(» - 1)(1 - /c) 2 (« + »c) V ' 



+ 



/c 4 (« 2 - 1) K*(n 2 - 1)(» 2 + c 2 ) 



/ 



12g BELL SYSTEM TECHNICAL JOURNAL 

The form of the final expression has been chosen so as to be suited to the 
use we shall make of it. 
Another double integral which appears in our work is 

,+« 
I(k, y) = / dv a(v ) 

(A2-10) 

fee 
dv a(v) exp [ — ik(v 4- v ) — y\v — vo\] 
.CO 

where a(v) is an even function of v and is such that all of the integrals 
encountered converge. We begin our transformation by dividing the 
interval of integration (- 00 , 00 ) for v into (- «> , 0) and (0, 00 ). Making 
the change of variable v = -v' , v = -v' in the first interval, dropping 
the primes and using a{— v) = a(v) leads to 

/(*, 7) = 2 f dv a(v ) f dv a(»)^ lf -^' cos k(v + v ) (A2-11) 

JO J— CO 

We now split the interval of integration of v in (A2-11) into the intervals 
(_ oo j 0), (0, v ), (v , °° ). In (- 00 , 0) we change the variable from v to 
-v', drop the prime, and use a{-v) = a(y). By paying attention to the 
sign of v — v we may remove the absolute value sign. By changing the 
order of integration in the double integral arising from the third interval 
(in which < v < 00 , » < < 00 ) we may show that it is equal to the 
double integral arising from the second interval. Thus 

j(k, 7) = 2 [ dv a(v ) f dv a(v)e~ yv ~ yv ° cos k(v - v) 
Jo Jo 

dv a(v ) / dv a(i)e- yV0+yv cos k(v + v) 
Jo 

When a(v), 7 and k are real we may write (A2-12) as 

/( fe >7 ) « 2 \f 'dva(v)e- yv - ikv2 
I Jo 

+ 4 Real f dv a(v )e-" l " +ikv ° f * dv a(v) e yv+ikv 
Jo J o 

and when 7 = ik we have 

I(k,ik) = 2 [ dv a(v ) [ dva(v)e- 2ikv 
Jo Jo 

+ 2 f dv a(v ) f° dv a(i)[e 2ikv + e~ 2ikv % 
Jo Jo 



(A2-12) 



(A2-13) 



(A2-14) 
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APPENDIX III 

Integral Equation When Guides Entering and Leaving Irregularity 
Are of Different Sizes 

Here we shall indicate how the integral equation method may be extended 
to cover the case mentioned in the above title. It is supposed that only the 
dominant mode is propagated freely in both guides. 

E in Plane of Irregularity 

Let the notation for the guide carrying the incident wave be the same as 
for the .E-corner. b denotes the narrow dimension of the guide and the 
quantities k and y m are given by (2-1). Both guides have the same wide 
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Fig. 3 

dimension a. The narrow dimension of the guide shown on the right of 
Fig. 3 is bi . We introduce the new quantity 

kx = [(2^/Xo) 2 - (Ai/a) 2 ] m (A3-1) 

to correspond to k. Since, by assumption, only the dominant mode is 
freely propagated in both guides both k and ki are real positive quantities 
less than unity. 

Let z = f(w) carry the system of Fig. 3 into a straight guide of width t 
in the w = v + id plane (see Fig. 2), and let g(v, 6) be denned by 

l + g(v,0)= \f'(w)\*. 

The behavior of g(v, 8) at infinity is shown by the table 

v dz/dw g(v, 6) 

— oo b/ir 

+ °° bi/ir k\W — 1 

where bi/b = k\/k has been used. It is convenient to introduce the ap- 
proximation g~(v) to g(v, 0). g~(v) may be chosen at our convenience subject 
only to the conditions that it be differentiable, |(— °o) = 0, and £(<») = 
klk~ 2 - 1. 
When we define G by equation (3-3) so that, as before, it is the Green's 
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function corresponding to a guide of width b, we may use equation (3-4) 
to derive the new integral equation 

Q(vo, 0o) = e- ikvo + £ f dvf dd 

LIT J- oo JQ 

• [g(v, 6)Q(v, 6) - £W E e- iklV }G(v , d ; v, B) (A3 " 2) 
+ T E F(vo) 
in which 

FW - r*'"»|W/g(») - e- ikeo N-(v ) - «**JV + (*) 

N-(v ) = I+kfa - k)- 1 [^ g'(v)e- i{ki - k)v dv (A33) 

N+(v ) = 2- I ft(*i + A)" 1 f tf(t)«" (fc,+,,) ' dv 

•>v 

Here g'(i>) denotes dg(v)/dv. Equation (A3-2) and 

Limit Q(v,6) = rj«-* lf (A3-4 } 

are to be solved for the unknown function Q(v, 0) and the unknown quantity 
T E • The method of successive approximations may be used in somewhat 
the same fashion as in the simpler case but we shall not give a general 
discussion. 
The first approximations are found to be 

A" = l/#-(oo), R? = -#+(-«>)/#-( oo) (A3-5) 

where the N's may be obtained by setting v = ±«i in equations (A3-3). 

One of the simplest choices for g(v) is to let it be zero for negative values 

of v and to have the value £(«>) = klk- 2 - 1 for positive values of v. Then 

T E l) = 2*(*i + /fe)" 1 , itf } = (* - *i) (* + ^i)- 1 (A3-6) 

These are quite similar to the corresponding expressions for a transmission 
line which have been used extensively in wave guide work. 

In working with these formulas, when k is small, it is sometimes convenient 
to use the result 

P dv f ' d6g{v,d) = r 2 r J P dv f d0\f (w)f - (tj - vj* (A3-7) 

J Vl ->0 J »l J o 

where the evaluation of the double integral on the right is made easier by 
the fact that it represents the area in the original guide (in the (x, y) plane) 
enclosed by the lines corresponding to v = Vi and v = v 2 . v 2 and vi are 
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chosen to be moderately large positive and negative numbers, respectively. 
It turns out that, when k\ and k are very small, this is related to the "excess 
capacity" localized at the irregularity whose effect must be added to that 
of the mismatch, indicated by (A3-6). 

When the entering and leaving guides are of the same size it is still possible 
to use the formulas of this appendix. N~(vo) may be replaced by an expres- 
sion which now has for its limiting value 

N-(oo) = 1 + i(k/2) [ g(v) dv (A3-8) 

H in Plane of Irregularity 

Let the figure corresponding to the irregularity be Fig. 3 with b and bi 
replaced by a and a\ , respectively. In addition to the quantities c and k 
defined by equations (5-2) we define 

K , = 2 fll /Xo , c, = {k\ - 1) 1/2 (A3-9) 

where we assume k and k\ to lie between 1 and 2. At v = — oo P(v, 6) 
still consists of the unit incident wave plus the reflected wave given by the 
first of equations (5-4) and g(v, 6) is still zero. However, now, at v = w , 

P(y, 0) = T,,e~ icxv sin 6 

£(oo) = kIk~ 2 - 1 = k~\c\ - c 2 ) (A3-10) 

The integral equation for P(v, 6) and T H is 

P(vo , 0o) = e~ icV0 sin 6o 4- ~ [ dv l 

Air J- » Jo 

•dO[g(T, 6)P(v, 6) - g{v)T H e- iciV sin d}G(v , 9, ; v, d) (A341) 

+ T a sin OqF h (vo) 
in which 

F H (v ) = <r'' Cl " g(r ) /£(«>) - e- ie ">M-(v ) - e icr ° M + (v ) 

•«'e 

(A3-12) 



M~(vo) = k 2 (2c)- 1 { Ci - c)~ l I* g'{i)e- i{c *- c) ° dv 

J— oo 



M + (v ) = « 2 (2 C )- 1 ( C + Ci r' f° g'(i)e-' Hc * +c) " 



</i' 



First approximations are 

Tg* = l/M-(oo), Rp = -M+(-oo)/M-(oo) (A3-13) 
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which, when we choose g\v) to be zero for v < and kik~~ — 1 for v > 0, 
become 

Iff = 2e( Cl + c)~\ R? = (c - o)fo + <0" 1 (A3-14) 

which again agrees with results obtained from transmission line considera- 
tions. When the entering and leaving guides are the same size we may use 

Jtf-(oo) = 1 + Ik (2c) _1 [ i(v) dv (A3-15) 

•Loo 

It seems difficult to give any general rules for the choice of g(v). Since, 
for R H and T a , the factor sin reduces the effect of the singularities on the 
walls of the transformed guide, the choice &(v) = g(v, sr/2) suggests itself. 
The factor sin 6 is not present in the formulas for R B and T B and regions 
near the walls are more important. In this case the selection 



i(v) = 7t" 1 f g (v,e) 

Jo 



dd 



may be useful, especially since it allows us to use the result (A3-7) when k 
and k\ become small. 

APPENDIX IV 
Variational Expressions for Reflection Coefficients 

The reflection coefficients are proportional to the stationary values of 
certain forms associated with the integral equations. In order to obtain 
these forms we proceed as follows. It is readily seen that the values of 
xi and # 2 which satisfy the symmetrical set of equations 

anxi + aux 2 = bi 

(A4-1) 

are the ones which make 

/ = dux* -(- 2auXiX2 + 022X-1 — 2biX\ — IbiXi (A4-2) 

stationary when xi and x 2 are given small arbitrary increments. This 
stationary value of J is 

J, = —blXl — &2*2 

If we take the integral equation to be the analogue of the set of linear 
equations, the reflection coefficient turns out to be proportional to 7, . 
In order to set down the actual expressions it is convenient to write r for 
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(v, ff) and dS for the element of area dvdd so that the integral equation 
(3-5) for Q(v, ff) may be written as 

Q(r ) = «-*" + ^ 2 (2tt)- 1 f g(r)Q(r)G(r , r) dS (A4-3) 

where the integration extends over the interior of the guide and G(r , r) 
denotes the Green's function (3-3). 

If the number of equations in the set (A4-1) were increased from two to a 
large number 2V, the set of x's would correspond, say, to the values of Q(r) 
or of g(r)Q(r), and the b's would correspond to the values of exp(— ikv ). 
In any event, we take the analogue of J to be 

Jb - / g{r)Q{r)[Q(r) - le'*"] dS 

(A4-4) 
- #(2t)- j ff g(r)Q(r)g(r )Q(r )G(r , r) dS dS 

where the subscript E indicates that we are dealing with an electric corner. 
It may be verified,* by giving Q(r) a small variation 5Q(r), that the function 
Q(r) which makes J E stationary is the one which satisfies the integral 
equation (A4-3). Furthermore, when we assume Q(r) to satisfy the integral 
equation, the expression for J E reduces to an integral which is proportional 
to the integral (3-6) for the reflection coefficient R B • More precisely, 
R E is given by 

ik 
R B = — [Stationary value of J E ] (A4-5) 

It follows that if, by some means, we have obtained a fairly good approxi- 
mation to Q, we may obtain a better approximation to R B by computing 
J E and using the formula 

R E = ik{2ir)-V E 

When we use the first approximation exp(— ikv) for Q to compute J B it 
turns out that the above formula gives the third approximation, R s \ to 
the reflection coefficient. 

The magnetic corner may be treated in much the same way. The 
integral equation (5-6) for P(v, ff) becomes, in the notation of this appendix, 

P(r ) = e - icv ° sin e + K 2 (27r) _1 f g(r)P(r)G(r , r) dS (A4-6) 

in which the v in dS = dvdd is integrated from — °o to -f- oo and 6 from to r, 

* See Courant and Hilbert, Methoden der Mathematischen Physik, Julius Springer, 
Berlin (1931), page 176, where a similar problem is treated. 
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as before, and G(r , r) now denotes the Green's function (5-5). We define 
Jnby 

Jh = [ g(r)P(r)[P(r) - 2e~ icv sin 0}dS 

(A4-7) 

~«W* ff i(r)P(r)g(r )P(ro)G(r Q , r) dS Q dS. 

J H is stationary with respect to small variations in P(r) when P(r) 
satisfies the integral equation (A4-6). Furthermore, from the integral 
(5-7) for R H , 

R l{ = u 2 (wc)- 1 [Stationary value of J H ] (A4-8) 

which may be used in the same way as equation (A4-5) for R E ■ 

J. Schwinger has used variational methods with considerable success to 
deal with obstacles in wave guides.* However, his variational equations 
differ somewhat from those given here. Some light on the relation between 
Schwinger's equations and the present one may be obtained by returning 
to the simple algebraic equations (A4-1) and (A4-2). A rough analogue 
of the expression required to be stationary in Schwinger's theory is 

(auxl + 2a 12 s 1 z 2 + a i2 xl)/(biXi + ta) 2 (A4-9) 

The essential point here is that the stationary value of the expression 
corresponding to (A4-9) gives the value of an impedance or combination 
of impedances appearing in some equivalent circuit. Expression (A4-9) 
may be obtained by expressing J, defined by (A4-2), as a function of xi 
and y = x 2 /xi . J is still to be made stationary but now it is a function of 
xi and y. Solving dJ/dxi = for x% and setting this value of xi in / gives 
the following function of y 

-(bi + M 2 fan + 2a ^y + *2*y*T~ 1 > 

which is the stationary value of / with respect to variations in Xi when y is 
held constant. This function is still required to be stationary with respect 
to y. The same is true of its reciprocal which becomes (A4-9) when both 
numerator and denominator are multiplied by Xi and the definition of y 
used. When (A4-1) is replaced by a larger number of equations similar 
considerations lead to a generalized form of (A4-9) . The expression required 
to be stationary by Schwinger is obtained when the sums in the general- 
ized form are replaced by integrals. 

* An account of the method together with applications is given in "Notes on Lectures 
by Julian Schwinger: Discontinuities in Waveguides" by David S. Saxon. An account 
is also given by John W. Miles. 11 
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